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SUMMARY 


A  simplex  computation  for  an  arc-chain  formulation  of 
the  maximal  multi-commodity  network  flow  problem  is  proposed. 
Since  the  number  of  variables  in  this  formulation  is  too 
large  to  be  dealt  v;lth  explicitly,  the  computation  treats 
non-baslc  variables  Implicif’y  by  replacing  the  usual  method 
of  determining  a  vector  to  enter  the  basis  with  several 
applications  of  a  combinatorial  algorithm  for  finding  a 
shortest  chain  Joining  a  pair  of  points  in  a  network. 
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A  SUOOESTED  COMPUTATION  FOR 
MAXIMAL  MULTI-COMMODITY  NETWORK  PLOWS 


Introduction 

A  problem  of  some  Importance  In  applications  of  linear 
programming  Is  the  determination  of  maximal  multi-commodity 
flows  In  networks.  For  example,  some  of  the  linear  programming 
prooleras  which  have  been  proposed  recently  by  Juncosa  and 
Kalaba  In  their  studies  of  communication  networks  [3]  can  be 
cast  In  this  form.  Straightforward  application  of  the  simplex 
method  to  such  problems  Is  usually  not  feasible,  since  even 
small  networks  may  generate  linear  programs  which  are  too  large 
for  present  machine  capacity.  What  Is  needed  are  specialized 
computing  schemes  that  take  advantage  of  the  structure  of  such 
problems.  For  the  single  commodity  case,  various  easy  compu¬ 
tations  are  known  [l,3#43,  but  the  multi-commodity  problem  has 
remained  relatively  unexplored. 

Consideration  of  simple  examples  makes  It  appear  that  the 
multi-commodity  flow  problem  Is  considerably  more  complex  than 
the  single  commodity  one.  Certainly  the  nice  romblnatorlal 
features  of  the  single  commodity  case  are  loot  In  the  generali¬ 
zation  -  simplex  bases  (for  any  formulation  of  the  problem 
known  to  us)  are  not  triangular,  hence  addition  and  subtraction 
do  not  suffice  to  solve  such  problems  by  the  simplex  method, 
the  max  flow  min  cut  theorem,  true  for  single  commodity  net¬ 
works,  Is  false  [3/  simple-minded  modification  of  the 

labeling  process  [V]  seems  to  work. 
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The  purpose  of  this  note  Is  to  suggest  a  computation  which 
makes  some  use  of  the  structure  of  one  formulation  of  the  multi¬ 
commodity  problem  within  the  framework  of  a  simplex  computation. 
Por  this  particular  formulation,  the  matrix  of  the  linear  pro¬ 
gram  Is  the  Incidence  matrix  of  arcs  vs.  all  chains  Joining 
sources  and  sinks  for  the  various  commodities,  and  thus  the 
number  of  variables  la  too  large  to  be  dealt  with  explicitly. 

The  suggested  computation  treats  non-baslc  variables  Implicitly 
by  replacing  the  "pricing" operation  of  the  simplex  method  (l.e. 
the  determination  of  a  vector  to  enter  the  basis)  with  several 
applications  of  a  combinatorial  algorithm  for  finding  a  shortest 
chain  Joining  a  pair  of  points  In  i  network. 


1 .  Arc-chain  Formulation 

Let  A^,  • • • #  be  a  list  of  the  arcs  of  the  network, 

Ci*  ...,  On  a  list  of  all  chains  that  Join,  for  the  various 
commodities,  all  the  sources  for  a  commodity  with  all  sinks  for 
the  same  commodity,  and  let  A  -  (a^^)  be  the  m  x  n  Incidence 
matrix  of  arcs  vs.  commodity  chains: 


{1  If  Cg  contains  A^, 
0  otherwise. 


Thus,  for  example.  If  the  network  Is  that  of  Pig.  1,  with  sources 
P^,  sink  P^  for  one  commodity,  and  source  P^,  sink  P^ 

for  a  second  commodity,  the  matrix  A  Is  as  shown  In  Pig.  2. 
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Plg.  1 


Cl  Cg  Cg  Cq  C^q 


Commodity  1 


^11  ^12  ^13  ^14  ^15 
1  1 

1  1 

1  1 

1 

1  1 
1  1 
Commodity  2 


Pig.  2 


If  we  let  X,,  8-1,  n,  denote  the  amount  of 

commodity  flow  along  C^,  and  the  flow  capacity  of  A^, 

then  the  multi— commodity  maximal  flow  problem  Is  represented 
by  the  linear  program: 
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n 

(2)  maximize  5 

3-1  ® 

to  the  constraints 
n 

g^^rs^s  ^  *n+r  “  ^r’ 

Xi*  ...,  >  0. 

The  assumption  In  (2)  that  commodities  are  valued  equally  Is 
not  essential  to  the  method  we  propose,  as  will  be  clear  from 
our  discussion  In  the  following  section.  Another  thing  we  wish 
to  point  out  Is  that  It  Is  Immaterial  whether  the  problem 
Involves  directed  or  undirected  arcs.  Thus,  for  example,  If 
there  are  "one-way  streets,"  or  If,  In  a  communication  network, 
say,  It  Is  desired  to  place  an  upper  bound  on  the  number  of 
messages  that  can  be  transmitted  from  and  P^,  and  an 

upper  bound  on  the  messages  that  can  be  sent  from  P^  to  P^, 
one  considers  two  arcs,  one  from  P^  to  Pj ,  the  other  from 
Pj  ^1'  directed  chains  from  sources  to  sinks. 

Since  the  number  of  chains  Is  usually  very  large  In 
practical  applications,  the  arc-chain  formulation  of  the 
problem  might  seem  to  be  Impossible  to  deal  with  computationally. 
Indeed,  the  enumeration  of  all  chains  from  commodity  sources  to 
sinks  In  a  network  of  moderate  size  would  be  a  lengthy  task,  to 
say  the  least.  Fortunately,  there  Is  no  need  to  write  down  the 
entire  matrix  A,  since  the  selection  of  a  variable  entering 
the  basic  set  at  ajiy  stage  of  the  simplex  computation  (or  the 


subject 

(5) 
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recognition  that  a  basis  is  optimal)  can  be  accomplished 
without  explicit  knowledge  of  the  non-baslc  column  vectors  of 
A.  All  we  need  Is  the  basis  B  ■  (or  Its  Inverse),  a 

square  submatrix  whose  order  Is  the  number  m  of  arcs  In  the 
network,  to  compute  the  simplex  multipliers  a^Cr  ■  l,...,m) 
satisfying,  for  J  •  s^,  ...,  s^, 

m  fl  If  J  <  n 

Z,Vrj  - 

r^  r  rj  [o  If  J  >  n. 

We  can  then  find  a  vector  to  bring  Into  the  basis  (or  prove 
that  the  current  basis  la  optimal)  by  the  method  of  the  next 
section.  Once  such  a  vector  has  been  found,  determination  of 
the  vector  leaving  the  basis  Is  accomplished  In  the  usual  way. 

2.  A  Shortest  Chain  Algorithm 

Suppose  we  have  computed  the  In  (4)  corresponding  to 

a  particular  basis  B.  If  some  Is  negative,  then  the 

variable  x^ may  be  Introduced  Into  the  basic  set  with 
n+r 

possibly  an  Increase  In  the  form  (2),  that  Is,  the  unit  vector 
having  1  In  the  r-th  position,  zeros  elsewhere,  can  be  brought 
into  the  basis.  (It  may  be  that  this  vector  also  represents  a 
one-arc  chain  for  some  commodity;  In  this  case,  a  bigger  Increase 
In  (2)  might  result,  of  course,  by  taking  the  latter  Interpre¬ 
tation.  ) 

Assume,  therefore,  that  a  stage  has  been  reached  In  the 
computation  where  all  are  non-negative.  In  this  case,  the 

algorithm  described  below,  which  makes  no  use  of  the  full 
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incidence  matrix  A,  can  be  used  either  to  locate  a  column 
vector  of  A  (l.e.  a  commodity  chain  In  the  network),  that  may 
be  brought  Into  the  basis,  or  to  prove  that  the  current  basis 
Is  optimal. 

Let  us  Interpret  the  as  lengths  of  the  arcs.  We  wish 

to  find  a  chain  C„,  If  one  exists,  whose  length 
ro 

^  a  a 

A 

Is  less  than  one,  the  coefficient  of  x.  In  (2).  Thus,  It 
suffices  to  locate,  for  each  commodity,  a  shortest  chain  from 
the  commodity  sources  to  Its  sinks.  If  each  of  the  chains  thus 
selected  has  length  at  least  one,  the  basis  Is  optimal.  Other¬ 
wise,  a  column  vector  of  A  corresponding  to  one  of  these 
chains  may  be  Introduced  into  the  basis. 

The  problem  of  locating  a  shortest  chain  from  one  set  of 
nodes  to  another  set  of  nodes  In  a  network  can  be  reduced  to  a 
standard  transshipment  problem  ,  and  may  consequently  be 
solved  In  various  simple  ways;  see  [6],  for  example. 

The  algorithm  we  describe  Is  that  of  [2j  .  (In  [2],  the  problem 
Is  considered  to  be  that  of  finding  a  shortest  chain  from  one 
node  to  another;  to  reduce  our  problem  to  this  one,  simply  Join 
each  node  of  the  first  set  to  a  new  node  by  an  arc  of  length 
zero,  and  similarly  for  the  other  net.  We  shall  give  a  des¬ 
cription  which  does  not  Involve  this  device  explicitly,  however.) 

Let  the  set  of  sources  for  one  commodity  be  S,  the  sinks 
for  the  commodity  T,  and  suppose  the  nodes  of  the  network  are 
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P^.  Let  denote  the  length  of  the  arc  joining 

and  Pj,  l.e.  If  the  arc  joining  P^  ,  and  Pj  corres¬ 


ponds  to  the  simplex  multiplier  set 


(If  arcs 


are  directed,  then  vie  let  denote  the  multiplier  corres¬ 

ponding  to  the  arc  from  P^^  to  Pj,  hence  In  this  case  we  may 
have  ^  whereas  in  the  undirected  case,  ■  ^jl*^ 

Initially  assign  to  each  node  P^  a  number  as  follows: 


^1  “ 


,0  for  P^  S 
00  otherwise. 


Now  scan  the  network  for  an  arc  S'lch  that 


^  <  'j- 


Replace  by  If  such  an  arc  Is  found.  Continue 

this  process.  Eventually  no  such  arcs  can  be  found;  then  the 
number  represents  the  length  of  a  shortest  chain  from  S 
to  P^,  for  all  1.  In  particular,  the  smallest  for 

Pj^  ^  T,  la  the  length  of  a  shortest  chain  from  S  to  T.  Let 
be  the  smallest  ouch.  To  find  a  chain  from  S  to  T  of 
length  TT^,  look  for  an  arc  P^Pj^  such  that  Tj 
then  search  for  an  arc  P^Pj  such  that  ■  Wj,  and  so 

on.  Eventually  a  node  of  S  Is  reached,  and  the  desired  chain 
has  been  traced  out  (in  reverse). 

If  In  the  process  of  locating  shortest  chains  from  commodity 
sources  to  sinks,  for  the  various  commodities,  one  Is  found  of 
length  less  than  one,  we  recommend  that  the  corresponding  column 
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vector  of  A  be  Introduced  Into  the  basis  immediately,  rather 
than  repeating  the  shortest  chain  algorithm  a  n\imber  of  times 
in  order  to  use  the  usual  criterion  for  selection  of  a  vector 
to  enter  the  basis. 

We  point  out  that  the  reason  for  getting  rid  of  negative 
multipliers  before  using  the  shortest  chain  algorithm  is 

that  the  algorithm  may  not  work  if  arcs  have  negative  lengths. 

To  start  the  simplex  computation,  one  can  of  course  begin 
with  the  basic  variables  •••»  ^n+r*  corresponding  to 

the  zero  flow. 

3«  Concluding  Remarks 

Except  for  hand  computation  of  a  few  small  problems,  we 
have  no  computational  experience  with  the  proposed  method. 

Whether  the  method  is  practicable  for  a  problem  Involving,  say, 

50  nodes,  100  arcs,  and  20  commodity  source-sink  sets 
Si,Ti,,..,  820**^20'  ^  question  which  can  be  settled  only 

by  experimentation.  It  would  certainly  be  more  practicable  in 
this  case  tl  an  straightforward  application  of  the  simplex 
method  to  a  node-arc  formulation  of  the  problem,  since  in  the 
latter  formulation  there  would  be  roughly  1100  equations  in 
2100  variables,  and  hence  the  basis  matrices  would  be  much  too 
large,  whereas  in  the  suggested  method,  the  basis  matrices 
would  be  100  x  100,  and  at  most  20  applications  of  the  shortest 
chain  algorithm  would  be  necessary  on  each  simplex  iteration.  How 
many  simplex  iterations  might  be  required  is  another  matter, 
though.  The  incidence  matrix  A  for  such  a  problem  could  have 
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nany  thousands  of  columns.  On  the  other  hand,  there  would 
probably  be  many  column  vectors  of  A  dominated  by  others.  In 
the  sense  that,  for  a  given  commodity  (or  for  different  com¬ 
modities  In  the  equal  value  case).  If  one  chain  C  Is  a  subset 
of  another  chain  C,  then  C*  can  be  Ignored.  (For  Instance, 
the  chain  of  Pig.  2  dominates  C0  and  Cg.)  The  shortest 
chain  method  takes  care  of  such  dominances  automatically. 

A  more  serious  consideration  Is  how  to  handle  the  case  of 
limited  supplies  of  commodities  In  such  a  problem.  For  example, 
suppose  that  In  the  two  commodity  maximal  flow  problem  corres¬ 
ponding  to  the  matrix  of  Fig.  2,  there  Is  an  amount  a^  of 
commodity  1  at  an  amount  a^  of  commodity  1  at  &nd 

an  amount  a|^  of  commodity  2  at  P^^.  We  can  reduce  this  to  a 
problem  of  the  same  type  as  before  by  Introducing  three  new 
directed  arcs  and  nodes  as  follows:  A^  from  P|  to  P^  with 
capacity  a^^,  A^  from  P^  to  P^  with  capacity  a^,  and 
AJ  from  PjJ  to  P^  with  capacity  a^.  We  then  take 
as  sources  for  commodity  1,  and  PjJ  as  the  source  for  commodity 

2.  However,  In  the  hypothesized  large  network  with  20  commodities, 

20 

the  number  of  such  new  arcs  would  be  2  »  where  n.  Is  the 

1-1  ^  ^ 

number  of  nodes  In  S^,  and  since  each  new  arc  Increases  the 
size  of  basis  matrices  by  one,  this  might  take  the  problem  out 
of  range  of  present  computing  machines. 
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